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Abstract. Polynomials in differentiation operators are considered. Joint 
covariance with respect to Darboux transformations of a pair of such poly- 
nomials (Lax pair) as a function of one-field is studied. Methodically, the 
transforms of the coefficients are equalized to Frechet differential (first term 
of the Taylor series on prolonged space) to establish the operator forms. In 
the commutative (Abelian) case, as it was recently proved for the KP-KdV 
Lax operators, it results in binary Bell (Faa de Bruno) differential polyno- 
mials having natural bilinear (Hirota) representation. Now next example 
of generalized Boussinesq equation with variable coefficients is studied, the 
dressing chain equations for the pair are derived. For a pair of generalized 
Zakharov-Shabat problems a set of integrable (non-commutative) poten- 
tials and hence nonlinear equations are constructed altogether with explicit 
dressing formulas. Some non-Abelian special functions are introduced. 



1. Introduction. 

Investigations of general Darboux transformation (DT) theory in the case 
of differential operators 

n 

L = £a fc 3 fc (1) 

with non-commutative coefficients was launched by papers of Matveev [1] . 
The proof of a general covariance of the equation 

ipt = Li> (2) 

with respect to the classic DT (the shorthands tp' = dtp = ip x are used 
through the paper) 

m=ip'-^ (3) 
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incorporates the auxiliary relation 

a t = dr + [r,a], r = J2o a n B n (a), (4) 

where B n are differential Bell (Faa de Bruno [3]) polynomials [2]. The re- 
lation (4) generalizes so-called Miura map and became the identity when 
a = <p'(f)~ l , cj) is a solution of the equation (2). 

Such operators (1) are used in the Lax representation constructions for 
nonlinear problems. It opens the way to produce wide classes of solutions 
of the nonlinear problem. Examples of discrete, non-Abelian and non-local 
equations, integrable by DT was considered in [4], [5]. Some of them were 
reviewed and developed in the book [6] and intensely used nowadays [7]. 
The approach was recently generalized for a wide class of polynomials of 
automorphism on a differential ring [8]. 

A study of jointly covariant combinations introduces extra problems 
of the appropriate choice of potentials on which the polynomial coefficients 
depend [10]. This problem was recently discussed in [11], where a method of 
the conditions account was developed. Covariant combinations of (general- 
ized) derivatives and potentials may be hence classified for linear problems. 
In two words, having the general statement about covariant form of a linear 
polynomial differential operator that determines transformation formulas 
for coefficients (Darboux theorem and its Matveev's generalizations), the 
consistency between two such formulas yields the special constraints. For 
example, the second order scalar differential operator has the only place for 
a potential and the covariance generate the classic Darboux transformation 
for it. 

In scalar case such one-potential constructions have been studied in [12] 
and developed for higher KdV and KP equations [13]. It was found that 
the result is conveniently written via such combinations of differentiation 
operator and exponential functions of the potential as Binary Bell Polyno- 
mials (BBP) [14]. The principle is reproduced and developed in the Sec. 
2.1 of this paper to give more explanations. 

The whole construction in general (non-Abelian) case is more compli- 
cated, but much more rich and promising. The theory could contain two 
ingredients. 

i) The first one would be non-Abelian Hirota construction in the terms of 
the mentioned binary Bell polynomials. On the level of general formulation 
some obstacles appears, e.g. an extension of addition formulas [14]) to the 
non-Abelian case . 

ii) The second way relates to some generalized polynomials that could 
be produced as covariant combinations of operators with a faith that ob- 
servations from Abelian theory could be generalized. Namely the case we 
would discuss in this paper. 
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Even the minimal (first order in the <9-operator) examples of the ZS 
problems with operator coefficients contains many interesting integrable 
models. It is seen already from the point of view of symmetry classification 
[15]. So, the link to DT covariance approach allows to hope for a realization 
of the main purpose - construction of covariant functions, their classification 
and use in the soliton equations theory. 

We would begin from the example, using notations from quantum me- 
chanics to emphasize the non-Abelian nature of the consideration. The 
operators p and H could play the roles of density matrices and Hamilto- 
nians, respectively, but one also can think of them as just some operators 
without any particular quantum mechanical connotations. The approach 
establishes the covariance with respect to DT of rather general Lax system 
for the equation 

-ip t = [H,h(p)], 

where h(p) - analytic function, in some sense -" Abelian" , i.e. the function to 
be defined by Taylor series [19]. More exactly it is shown that the following 
statement takes place 

Theorem. Assume (\\ and {ip\ are solutions of the following (direct) 
equations 

-i(Xt\ = ~(x\Hp), 
v 

and \(p) stands for the conjugate pair. Here p, H are operators left- acting 
on a "bra" vectors (ip\ associated with an element of a Hilbert space. The 
transforms pi, h(p)i are defined by 

(^1 = ^1(1 + ^^), (5) 

Pl = TpT-\h 1 (p) = Th(p)T-\ T=(l + ^— -P). (6) 

where P = \ip)(x\ l f)(x\- Then the pairs are covariant: 

z\{i>\\ = (4>i\(pi ~ A-ff), = -(ViIMp)- 

complex numbers X, z\ are independent oft [19]. 

The cases f(p) = ip 3 and f(p) = ip~ x were considered in [22], see 
applications in [17]. A step to further generalizations for essentially non- 
Abelian functions, e.g. h(X) = XA + AX, [A, X] ^ 0, is studied in [16]. 
The case is the development of the matrix representation of the Euler top 
model [20]. This example of the theory is more close to the spirit of the 



4 



sec. 3.4, more achievements are demonstrated in [21], where abundant set 
of integrable equations is listed. The list is in a partial correspondence 
with [15], and give the usual for the DT technique link to solutions via 
the iteration procedures or dressing chains. One of the main results, we 
present in the sec. 3.3, is how the "true" non-Abelian functions appear in 
the context of the covariance conditions application. 

2. One- field Lax pair for Abelian case. 

2.1. COVARIANCE EQUATIONS 

First we would reproduce the "Abelian" scheme, generalizing the study of 
the example of the Boussinesq equation [11]. To start with the search we 
should fix the number of fields. Let us consider the third order operator (1) 
with coefficients bk, k = 0, 1, 2, 3, reserving a& for the second operator in a 
Lax pair. Suppose, both operators depend on the only potential function w. 
The problem we consider now may be formulated as follows: To find restric- 
tions on the coefficients bs(t),b 2 (x, t), b\ = b(w,t), bo = G(w,t) compatible 
with DT transformations rules of the potential function w induced by DT 
for 6j. The classic DT for the third order operator coefficients (Matveev 
generalization [1]) yields 

b 2 [l]=b 2 + b' 3 , (7) 
6i[l] = b 1 + b' 2 + 3b 3 a > , (8) 
b [1] = b + b[ + ab' 2 + 3b 3 (aa' + a"), (9) 

having in mind that the "elder" coefficient 63 does not transform. Note also, 
that b 3 = yields invariance of the coefficient b 2 . 

The general idea of DT form-invariance may be realized considering the 
coefficients transforms to be consistent with respect to the fixed transform 
of w. Generalizing the analysis of the third order operator transformation 
[11], one arrives at the equations for the functions b 2 (x,t),b(w,t),G(w). 
The covariance of the spectral equation 

b^xxx + b 2 (x, t)ip xx + b(w, t)ip x + G(w, t)ip = Xip (10) 

may be considered separately, that leads to the link between 6j only. We, 
however, study the problem of the (10) in the context of Lax representation 
for some nonlinear equation, hence the covariance of the second Lax equa- 
tion is taken into account from the very beginning. We name such principle 
as the "principle of joint covariance" [10]. The second (evolution) equation 
of the case is : 

ipt = a 2 (t)ip xx + ai(t)ip x + wip, (11) 
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with the operator in the r.h.s. having again the form of (1). We do not 
consider here a dependence of ctj , b% on x for the sake of brevity, leaving this 
interesting question to the next paper. 

If one consider the L and A operators of the form (1), specified in 
equations (10) and (11). as the Lax pair equations, the DT of w implied 
by the covariance of (11), should be compatible with DT formulas of both 
coefficients of (10) depending on the only variable w. 

a 2 [l] = a 2 = a(x,t), 

ai[l] = ai(x, t) + Da(x, t) 

a [l] = w [l] = w + a[ + 2a 2 a' + aa 2 (12) 

Next important relations being in fact the identities in the DT transfor- 
mation theory [2], are the particular cases of the generalized Miura map 
(4): 

a t = [a 2 (a 2 + a x ) + aicr + w] x (13) 
for the problem (11) and, for the (10) 

hi^ 3 + 3a x a + a xx ) + b 2 (a 2 + 
a x ) + b(w, t)a + G(w) = const; 

is a solution of both Lax equations. 

Suppose now that the coefficients of the operators are analytical func- 
tions of w together with its derivatives (or integrals) with respect to x (such 
functions are named functions on prolonged space [24]). For the coefficient 
bo = G(w, t) it means 

G = G(d^ 1 w, w, w x , d~ 1 w t , w t ,w tx , ...)■ (15) 

The covariance condition is obtained for the Frechet derivative (FD) of the 
function G on the prolonged space, or the first terms of multidimensional 
Taylor series for (15), read 

G(w + a[ + 2a 2 cr' + aa' 2 ) = G(w)+ , * 

G Wx (a' 1 + 2a 2 a' + aa' 2 )' + .... 1 ° J 

We shall show only the terms of further importance. 

Quite similar expansion arises for the coefficient b\ = b(w,t), with which 
we would start in the analogy with the expressions (8, 16). Equalizing the 
DT and the expansion one obtains the condition 

b 2 + 3b 3 a' = 6 w (ai + 2a 2 a' + aa' 2 ) + b w >(a[ + 2a 2 a + aa 2 )' .... (17) 
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This equation we name the (first) "joint covariance equation" that guaran- 
tee the consistency between transformations of the coefficients of the Lax 
pair ((10), 11). In the frame of our choice a 2 = 0, the equation simplifies 
and linear independence of the derivatives yields two constraints 

363 = 2b w ci2, ,.. q \ 

b' - b a 1 ^> 

or, solving the second and plugging into the first, results in 

b w = 3b 3 /2a 2 , hQ , 
b' 2 = 363a; /2a 2 . [ ' 

So, if one wants to save the form of the standard DT for the variable w 
(potential) the simple comparison of both transformation formulas gives 
for b(w) the following connection (with arbitrary function a(t)) 

b{w, t) = 3b 3 w/2a 2 + a(t). (20) 

Equalizing the expansion (16) with the transform of the feo = G(w,t) 
yields: 

V 1 + aV 2 + 3b 3 {o 2 /2 + <j , y = 
G Wx {a[ + 2a 2 a' + aa' 2 )' + G d -i Wt [a lt + 2d- 1 {a 2 a' t ) + d^iaa'^t] + ... 

(21) 

This second "joint covariance equation" also simplifies when a 2 = 0: 

3b 3 w'/2a 2 + ab' 2 + 36 3 (d~V t - w)'/2a 2 + 3b 3 a"/2 = G^a^ + 2a 2 a')'+ 
G d -i w [ai + 2a 2 a] + G d -i Wt [a lt + 2a 2 a t ] + 

(22) 

when (20) is accounted. Note, that the "Miura" (13) is used in the l.h.s. 
and linearizes the FD with respect to a. Therefore, the derivatives of the 
function G 

G Wx = 36 3 /4a 2 , 
G d -i wt = 36 3 /4a|, (23) 
Gq-^w = b' 2 /2a 2 , 

are accompanied by the constraint 

an + a 2 a" + axa\ = 0, (24) 

which have got the form of the Burgers equation after (19) account. Finally 
the integration of the relation (19) gives 

b 2 = 36 3 ai/2a 2 + 0(t) (25) 
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and the "lower" coefficient of the third order operator is expressed by 

G(w, t) = 3b 3 w x /2a 2 + 3b 3 a' 1 d- 1 w/2(a 2 ) 2 + 3b 3 d' l w t /2a 2 2 . (26) 

Statement 1 The expressions (11, 10, 20, 26) define the covariant Lax pair 
when the constraints ( 19, 24) are valid. 

Remark We cut the Frechet differential formulas on the level that is 
necessary for the minimal flows. The account of higher terms leads to the 
whole hierarchy [13]. 

2.2. COMPATIBILITY CONDITION. 

In the case a! 2 = by which we have restricted ourselves, the Lax system 
(10, 11) produces the following compatibility conditions: 

2a 2 b' 3 = 3b 3 a 2 , 
b 3t = 2a 2 b' 2 - 3b 3 a'{ 
b 2t = a 2 b' 2 ' + 2a 2 b' 1 + a x b' 2 - 3b 3 a'( - 2b 2 a[ - 3b 3 a' (27) 
ht = a 2 b'{ + ai&i - b 3 a'{' - b 2 a'{ - Mi - 36sog - 2b 2 a' + 2a 2 b' 
bot = aib' G + a 2 b f Q - Mo - b 2 a^ - b 3 a" 

In the particular case of a 2 = we extract at once from the first of the 
equalities (27) the constraint b' 3 = 0. The direct corollary of (25) is b 3 t = 0. 
In the rest of the equations the links (27,25) are taken into account. Hence 
(24) in the combination with the expression for b 2 t produce 

0t = -2(3a[. (28) 

The last two equations (choice of constants 63 = l,a 2 = —1) become 

aw + a t + 3a'{d- 1 w/2 + (2(3 - 3a 1 /2)w' + a'f + 3a 1 a'(/2 = 
39" 1 (w t + a lW ) t /A = (a - 3w/2)w' - w'" /4+ 
3a lWt /A + 3a 1 a'{d' 1 w/^ + 3aia>/4 - 3a>'/4+ [ ' 

(P + 3 ai /4)w". 

In the simplest case of constant coefficients (b' 2 = a[ = 0) one goes down to 

3b 3 (w t + a±w) t /4:a 2 = 
-\{3b 3 w/2a 2 + a)w' - b 3 w'"/4 + 3b 3 a 1 w t /Aa 2 2 + {(3- 3b 3 a 1 /Aa 2 )w"]' . 

(30) 

This equation reduces to the standard Boussinesq equation when (61 = 
ai = 0,63 = l,a 2 = -1) [6]. 

We would repeat that the results given in the Sec 2 are simplified to 
show more clear the algorithm of the covariant Lax pair derivation. More 
general study (a 2 / 0) will be published elsewhere. 
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2.3. SOLUTIONS. DRESSING CHAINS FOR THE BOUSSINESQ EQUATION 

The dressing formula for the zero seed potential (39) is standard and in- 
cludes the only seed solution cp, of the Lax equations with zero potential 
w. 

w s = a[ + 2a2<j' + aa' 2 = a[ + 2d2 \og xx (j)(x, t) (31) 

A next power tool to obtain solutions of nonlinear system is the dressing 
chain equation: solitonic, finite-gap and other important solutions were ob- 
tained for the KdV equations reducing such chain [23] . Going to the dressing 
chain, we use the scheme from [11]. We would restrict ourselves further to 
the case of a 2 = —1, ao = u, 63 = 1, b 2 = 0, 61 = b(u, t) = —3u/2 + a,G = 
—3u'/4 + 3d~ 1 ut to fit the notations from [11]. The general construction is 
quite similar. 

The Miura equations (13, 14) also simplifies: 

a t = -{<J 2 + <Jx)x + u x (32) 

for the problem (11) and 

a 3 + 3a x a + a xx + ba + G = const, (33) 

where b = 3u/2 + a, G = -3d~ l u t /A + 3u x /4. 

Namely the equations (32, 33) together with the n-fold iterated DT 
formula (5) 

u n+1 = u n - 2a' n (34) 

form the basis to produce the DT dressing chain equations. 
We express the iterated potential w n from (32) 

-CT nt + {a 2 n + a'J = u' n (35) 

and substitute it into the differentiated relation (34) to get the first dressing 
chain equation 

<Tn+l,t - CTnt = (o-n+1 + a n+lY ~ i a l ~ ( 36 ) 

Next chain equation is obtained when one plugs the potential from (35) to 
the iterated (33) 

al + 3a' n a n + a^ + (-3u n /2 + a)a n + -3u' n /A + 3d~ 1 u n t = c n . (37) 
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3. Non-Abelian case. Zakharov-Shabat (ZS) problem. 

3.1. JOINT COVARIANCE CONDITIONS FOR GENERAL ZS EQUATIONS 

Let us change notations for the first order (n=l) equation (1) with the 
coefficients from a non-Abelian differential ring A (for details of the math- 
ematical objects definitions see [2]) as follows 

ipt = (J + ud)ip, (38) 

where the operator J G A does not depend on x, y, t and the potential 
ao = u = u(x, y,t) € A is a function of all variables. The operator d = d/dx 
may be considered as a general differentiation as in [2]. The transformed 
potential 

u = u+[J,a], (39) 

where the a = (f) x (f)~ l : is defined by the same formula as before, but the 
order of elements is important. The covariance of the operator in (38) fol- 
lows from general transformations of the coefficients of a polynomial [6]. 
The coefficient J does not transform. 

Suppose the second operator of a Lax pair has the same form, but with 
different entries and derivatives. 

ip y = (Y + wd)^, (40) 

Y € A where the potential w = F(u) G A is a function of the potential of 
the first (38) equation. The principle of joint covariance [10] hence reads: 

w = w + [Y,a}= F(u+[J,a}), (41) 

with the direct corollary 

F(u) + [Y,a}=F(u+[J,a}). (42) 

So, the equation (42) defines the function F(u), we shall name this equation 
as joint covariance equation. In the case of Abelian algebra we used the 
Taylor series (generalized by use of a Frechet derivative) to determine the 
function. Now some more generalization is necessary. Let us make some 
general remarks. 

An operator-valued function F(u) of an operator u in a Banach space 
may be considered as a generalized Taylor series with coefficients that are 
expressed in terms of Frechet derivatives. The linear in u part of the series 
approximates (in a sense of the space norm) the function 



F(u) = F(Q) + F\0)u + 
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The representation is not unique and the similar expression 

F(u) = F(0) + uF'(0) + ... 

may be introduced (definitions are given in Appendix). Both expressions 
however are not Hermitian, hence not suitable for the majority of physical 
models. It means, that the class or such operator functions is too restrictive. 
To explain what we have in mind, let us consider examples. 

3.2. IMPORTANT EXAMPLE 

From a point of view of the physical modelling the following Hermitian 
approximation 

F(u) = F(0) + H + u + uH + u + = u, 

is preferable. Such models could be applied to quantum theories: introduc- 
tion of this approximation is similar to "phi in quadro" (Landau-Ginzburg) 
model [16]. Let us study, in which conditions the function 

w = F(u) =Hu + uH, (43) 

satisfy the joint covariance condition for the Lax pair (38, 40) By direct 
calculation in (42) one arrives at the equality 

[Y,a]=H[J,a] + [J,a}H. (44) 

The obvious choice for arbitrary a is Y = H 2 , J = H. 

The compatibility conditions for the pair of equations (38) and (40) 
yields 

u y - Hu t - u t H + [u, H]u + u[u, H] + H 2 u x + Hu x H + H 2 u x = (45) 
If the potential does not depend on t, it is reduced to the next equation: 
u y + [u 2 , H] + H 2 u x + Hu x H + H 2 u x = 0, (46) 
and x-independence yields the generalized Euler top equations 

u y + [u 2 ,H] =0, (47) 
which Lax pair (38,40) with Y = J 2 , J = H was found by Manakov [20]. 

3.3. CO VARIANT COMBINATIONS OF SYMMETRIC POLYNOMIALS 
The next natural example appears if one examine the link (44) . 

P 2 (H, u) = H 2 u + HuH + uH 2 



11 



The direct substitution in the covariance and compatibility equations leads 
to covariant constraint that turns to the identity, if Y = H 3 , J = H. 

It is easy to check more general connection Y = J n ,J = H connection 
that leads to the covariance of the function 

P n (H,u) = J2H n - p uHP. 

Such observation was exhibited in [16]. On the way of a further generaliza- 
tion let us consider 

f{H,u) = Hu + uH + S 2 u + SuS + uS 2 (48) 

Plugging (48) as F(u) = f(H,u) into (42), representing Y = AB + CDE 
yields 

A[B, a] + [A, a]B + CD[E, a] + C[D, a]E + [C, a]DE = H[J, a] + [J, a]H+ 
S 2 [J,a}+ S[J,a]S + [J,a]S 2 . 

The last expression turns to identity if A = B = J = H,C = aH, D = 
ocH, D = aH, S = 0H and [a, H] = 0, \j3, H] = with the link a 3 = (3 2 . 

Statement Darboux covariance define a class of homogeneous poly- 
nomials P n (H,u), symmetric with respect to cyclic permutations. A lin- 
ear combination of such polynomials J2n=i PnPn(H,u) with the coefficients 
commuting with u, H is also covariant, if the element Y = J2n=i d n H n+1 
and ai = 0i = 1, <+ 2 = P n+1 , n + 1. 

A proof could be made by induction that is based on homogeneity of 
the P n and linearity of the constraints with respect to u. The functions 
Fh{u) = X)o° a nPHn{u) satisfy the constraints if the series converges. 

4. Conclusion 

The main result of this paper is the covariant equation (42). See, also the 
example (44). A class of potentials from [21] contains polynomials P n (H,u) 
and give alternative expressions for it. The linear combinations, introduced 
here could better reproduce physical situation of interest. So, we used the 
compatibility condition to find the form of integrable equation and reduc- 
tion tracing the simplifications appearing for the subclasses of covariant 
potentials. While doing this we also check the invariance of the equation 
and heredity of the constraints. 

The work is also supported by KBN grant 5P03B 040 20. 
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derivative of an operator by other one is defined in a Banach space B Two 
specific features in the case of a operator-function F(u) <G Bu G B should 
be taken into account: a norm choice when a limiting procedure is made and 
the nonabelian character of expressions while the differential and difference 
introduced. 

Definition. Let a Banach space B have a structure of a differential 
ring. Let F be the operator from B to B' defined on the open set of B. The 
operator is named the left-differentiable in uq E B if there exist a linear 
restricted operator L(uq), acting also from B to B' with the property 

L(u + h)- L{uq) = L(u )h + a(u Q , h), \\h\\ -> 0, (49) 

where ||a(«o, /i)||/||/i|| — > 0. The operator L(uq) = F'(uq) is referred as 
the operator of the (strong) left derivative of the function F(u). The right 
derivative F'(uq) could be defined by the similar expression and conditions, 
if one changes Lh — > hL in the equality (49). 

The addition of the half of the right and left differentials 

(F'(u )h + hF'(u ))/2 (50) 

also approximates the difference L(no + h) — L(uq) in the sense of the FD 
definition. 



